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Abstract 

The mean square radius of the proton charge distribution was studied in the 
framework of the relativistic quasipotential quark model in assumption of the SU(6)- 
symmetry. It was shown that the proton charge radius is represented as function 
/(m/7)/M^ for point quarks (m is a quark mass, 7 is the scale of the coupling energy, 
M is the nucleon mass, / is undimentional function, which does not depended on 
M). It was shown, that in the ultra relativistic region the mean square radius of 
the bound system can have negative value. To describe simultaneously nucleon 
magnetic moments and proton radius in the oscillator model it is need to suppose 
that quark has the negative mean square radius (rf^) = — 1.875/m^. 

Introduction. 

In works [1, 2, 3] the formalizm of the quark model in the quasipotential approach 
was considered: the nucleon quark wave function was constructed in the impulse approx- 
imation [1], the calculations of the magnetic moments [2] and the ratio of the axial and 
vector coupling constants [3] were performed. It was obtained that model has property 
which is followed by the relativistic kinematics of quarks: nucleon magnetic moments (in 
units of nuclear magnetons) and ratio Ga/Gv depend on one undimensional parameter, 
which is constructed as ratio of the quark mass and the energy scale parameter, and don't 
depend on the ratio of the quark mass and nucleon mass. In this case the valent qaurk 
mass and energy of coupling are intrinsic parameter, which are related. Nucleon charge 
radius, considered in the this paper, has this property. 

1. Basic relations of the model 



X 
J3 



In the quasipotential quark model [1] the electromagnetic current in the rest frame of 
the initial nucleon has the form [2]: 

(KA'^A^ I J^(0) I OAjAt) = eJ^'/^-^(K,0)5A' A, = 



3e / dnp^dn^^^{MQ)Gf,^{MQ)5r>r, 



where 



^ Pi P2 P3 



Mo = Ep^ + + i?p3 , 



Eo'E. 



Pa 



P3 



r(1) 
r(2) 

r(3) 

.(3)A' As^, , , 

e> (Ps.PsJ'^Ai^Atg 



-'^A',, A,, l-'^Pi2'Pl J-'^A^,2A,2l-^Pi2> P2 , 

(P3Al3Ayjf (0)|P3A.,A,,) 



(2) 
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Here jjj^\0) is a current of third quark. For convenience the electron charge e is written 
as factor. Sum of 3-momenta in the wave function of the relative motion (p is equal to 

o o o ' 

zero: Pi + P2 + P3= 0, Pi + P2 + Ps = 0. The list of notations is given at the end of 
the section. The current normalization condition Jo(0, 0) = ejv (ejv is the nucleon charge 
in units of the electron charge e) leads to the WF normahzation condition [1], which is 
given by: 

J dnp.dQp, I (^m(Mo) 1 2 /Ep, = 1. (3) 
The Dirac form factors of the nucleon are defined by the expression: 

Ji'^^^(K, 0) = er^U'Hm^.HQ') + ^HQ>,uQu}U'^iO), (4) 



where k is the nucleon anomalous magnetic moment, cr^jy = ^{'Jt^'Ju—^u^fj), = —Q^ > 0, 
q — (P' — P), P' = (E^jK.), P = (M, 0). Sachs form factors is convenient for analysis 
of the experimental data: ^^^(Q^) = Fi{Q^) + F2iQ^)KQyAM\ GuiQ^) = Fi{Q^) + 
K,F2{Q'^). Let write the zero component of the third quark current like (4), but in form 
of Pauly cr— matrix, with zero quark anomalous magnetic moments: 

(P3, Pa) = ^MQh) + + ^) + {ap',){ap,)], (5) 

2^{E3 + m){E'3 + m) 

where 63 is charge operator of the third quark in units of the electron charge, /i is the 
quark Dirac form factor, depended on Q^g^ = — (^3 — Ps)^- Quark radius is defined by 
decomposition over Qf^y. fkiQf^)) ^ 1 - {rl^)Q%)l^- 

In work [1] the relativistic three-particle oscillator was proposed in the effective mass 
approximation. This approximation is defined by the condition: rn^J^ — (Ep^) « nfj^ , 



where m^^-^ = ^w? + (p^), (p^) = S d^p^d^p^\ipMVPk/ Ep.,, {EpJ = S d^p^d^p^\ipMVEpJ Ep^. 
The wave function (/? of relative motion in this approximation is written as: 

o m 3k'^ 

V'M ({Pfe}) = A^exp[--(Mo - 3m)] « A^exp(-- - ^) ^ (6) 

1 o2 o2 00 2m . 

^ TVexp -— Pi + P2 + P1P2 ^ ^ff ). 



where k = (tti — vr2)/2, k' = (tti + 7r2)/3 — 27r3/3; tt^ =Pfc ,/2m/ im + Eo ) is half- 

V Pfc 

momentum [4]; the energy has the form Eo = m + 7r^/2m. In the effective mass limit 

Pfc 

O I O 2 

T^k ~Pfc V 2m/ (m + m^ff), tti + 7r2 + tts 0, Eo m+ Pj. /(m + rrf^^). If (pi) is small, 
then m^^^ ~ m and effective mass approximation passes to nonrelativistic limit. 



Here we give a list of notations: dVlp^ = dpk/Ep^,, Ep^ = ypf. + m^, E^ = + M^, 
K and are nucleon and quark momenta, M and m are nucleon and quark masses, U 
and u are nucleon and quark spinors, {pkKk^t^} = Pi-^siAtiP2A<j2At2P3As3At3; x is SU(3) 
nucleon wave function. The total spins and isospins of the nucleon and constituent quarks 
are omitted and equal to 1/2; Aj, A^^^, are the third projections of the nucleon and quark 
spins and equal to +1/2. We assume a sum over the repeated indices of the third pro- 
jection of spin and isospin. The normalization condition of vector states and spinors of 
nucleon and quarks are given by: 

{K'J'X'jT'X'^ I KJXjTXt) = (27r)353(K' - K)5,. 
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M^^fe(p)M^^fc(p) = m, U^' {K)U^-' (K) = M/E^ (there is no sum over A,, and Aj). 



Pk= L-^Pk, where pi2 = Pi +P2; Pfe= (VPfe), ^^fe), where Pi2=Pi + P2; 



P12 



is the Lorentz boost to the rest frame of the nucleon and L ^ and Lz^^ are the Lorentz 

boosts to the rest frame of [1+2] quarks. Wigner rotation matrix has the form (see [1] 
and references there in): 



D^"'\LK\Pk) 



(Ep,+m)(E^ + M)-(aK)(ap,) 
2(Ep, + m){E^ + M){E^E^^ - Kp^ + mM) ' 



(7) 



4-momentum in the c.m.s. has the components: 



Kpfe), Pfe=Pfe 



K Kpfe 



(8) 



3. Analitycal analysis. 

By using the current (5), the expression for Ge{Q'^) is represented as: Ge{Q'^) 



2-E'k / + -E^K )>^o(K, 0). Electric nucleon radius, defined according to the decomposi- 
tion GEiQ"^) ~ ejv — {r%)Q'^/6, expresses as: 



(4) = -v^GE(g2)|K=o 



3e 



AT 



4M2 



-V^Jo(K,0)||K=o. 



(9) 



In the nonrelativistic limit we have from (9) and (1), denoting ejv = 3(x, e^^^x): 

{rlf ^CnJ rfpirfp2(-iVK) V(Pi, P2, P3)¥'(Pl, P2, P3). (10) 

Using the last expression of (6) with m'^^-^ ~ m, we have from (10): 



2\NR 



Ctv j (ipi(ip2|97(pi,P2,P3) 



(11) 



Now if we return to variables k and k', assuming pi = k + k'/2, p2 = — k + k'/2 and 
using expression for the relative motion wave function (6) in variables k and k', wc have 
well-known results [5]: {r\)^^ = where ortonormalization condition (3) is used 

also. Now we can rewrite general expression for radius in the relativistic case with the 
oscillator wave function (6): 



(4) 
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4M2 



m 



m , 
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-eiv-T(VK MoY-ET- + eiv-^VK' Mq "f^ + 



7 



E, 



P3 



7 



E, 



P3 



m 



+3^(Vk Mo)(VkG'o) - 3Vk'Go 
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|K=0 



eiv[^ + {r}) + {rh) + {rhi) + (r,V)]- (12) 



Four terms {rj u m jy) correspond to four expressions under integral in the expres- 
sion (12). After differentiation over K for K = we don't have any cr-matrix, that 
is why all expressions are proportional to the nucleon charge = 3(x, Cg^-'x)- So, in the 
quark model the neutron radius equal to zero, but neutron form factors deviate from zero 
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for unzero transfer momenta [6]. Further we consider only the proton radius. Expressions 
for (rf ) and (r|j) have the forms: 

^^i) - M-^^^J ^^^Pi^^^P^I^I E„, El w? 



M2 ^- ^ ^ ^ 



3{Ep^ + Ep^) + ^(-Ep3 - Ep^ - Ep^] 
%3 



T.Ep^ (14) 



Ep^m 



■"P3 

Third term {r'jjj) contains derivative from Go (2). In this expression B'^^^ does not depend 
on K. It was shown numerically that expression Vk-B*^^^|k=o is negledibly small in the 

wide regionof the parameters and Vk-B*^^-* |k=o is proportional to pi + p2 =0 in the 
effective mass limit. Taking into account that Pi + P2 = — Pa, it is possible to show that 
contribution of B^'^^ is equal to zero. So, we have only contribution from I/-E0' and from 

the quark current: 

Nonrelativistic expression for radius (eAr/7^) can be obtained from (13)-(15), assuming 
M pa 3m. 

Relativistic term (r|y) is represented as series of the negative powers [m/7]^ (including 
the zero power). As a result we have: 

9 4 2 2 4 

= TT^(—a-4 + —a-2 + ao + — 02 + — 04 + ■■■)■ (16) 

\ ^4 777,^ 777,4 y 

Therefore we don't have dimensional terms in the form 1/m^ for radius (16) (compare [5, 
7]). 

Dependence on the nucleon momentum K in (1) is defined by Lorentz-transformation, 
which contains 3-vclosity K/M (see (8) ). Therefore, the nucleon radius (9) for arbitrary 
wave function is defined as /(m/7)/Af^, where / is undimensional function.If quarks have 
structure, then / depends on the parameters of the quark structure. 

3. Numerical analysis. 

Figure a shows dependences of proton radius {r^) (curve 1), (r|) + {r}j) + {r'jjj) 
(curve 2) and (r|y) (curve 3) from m/7. Points for {r\) was obtained as follows: data 
for current and for electric form factor Ge{Q'^) were obtained for fixed m/7 and in the 
interval 0.001 < < 0.1 GeV^. Integral was calculated by VEGAS [8] with five iteration 
and 10^ trials in each. After that points of Ge{Q^) were fitted by function Ge{Q^) ~ 
bo — 6iQ^/6 + b2Q'^ + b^Q^ + b^Q^ for obtaining the radius {r^) = bi. Errors of the fit are 
shown on the figures. 

For calculation of (r|) + (r|j) + (r|jj) we can use (13)-(15), calculated by VEGAS [8]. 
Error of the calculation is small and shown on the figure. For calculation of (r|y) it was 
difficultly to obtain analitycal formulae, like (13)-(15). However, we can separate the 

o I 

contribution of (r^y) in the current (1). For this it is sufficiently to replace ^^(Mo) 
(P{Mq) and perform calculation as in case of (r|;). 

If m/7 — > 00 then (r|;) increases because of there are positive powers of m/7, as in 
nonrelativistic theory. If m/7 decreases then relativistic term with negative powers of 
m/7 dominates. 
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Prom figure a we can see that (rj) + (rjj) + {rjjj) has nonzero hmit for m/7 — > 0. 
Indeed, it is possible to represent integrals in (13), (14), (15) as infinite series of the 
decomposition with the negative powers of m/7. This decompositions have nonzero terms 
[m/7]~*^ in (13) and [m/7]~^ in (14) and (15), which define limits m/7 — > for (r|), {rjj) 
{rjjj). Nonzero limit of strong coupling for the proton radius is relativistic effect. In the 
nonrelativistic case size of system decreases to zero when the bound energy increases. 

Prom figure a we can see, that behaviour of the radius of the bound system is more 
complicated in the limit of strong coupling as compare to the behaviour of (rj) + {rjj) + 
(rjjj). In the region of small m/7 form factor GEiQ"^) > 1 for small Q^, radius of the 
bound system is negative and its module increases when m/7 — > 0. We can analyze 
what terms in the expression (1) correspond to this behaviour of the form factor in 
the ultrarelativistic region. From figure a we can see that the cause is the term {rjy). 
According to (1), (2) contribution of (rjy) consist of following parts: contribution of the 
Wigner rotation matrix {r]^) and contribution of the quark current (r|). Cross term in 
which first derivative Vk acts to Wigner rotation matrix, second derivative acts to the 
quark current, does not give any contribution (see explanation before (15)). Figure b 
shows the proton radius dependences on m/7: (r^) (curve 1), (r|) (curve 4), (r|,) (curve 
5). It is possible to conclude, that the positive peak and the negative proton radius for 
small m/7 is originated from the quark current. So, the reason for the negative radius of 
the bound system of the spinor particle is the dependence of the valence quark current 
on quark momenta, defined by expression (5). 

The quark radius contributes to the terms with negative powers of m/7 in the ex- 
pression (16). We can write down the quark radius contribution to the proton radius 
exactly. 

In work [2] for oscillator model (6) it was shown that simultaneous description of proton 
and neutron magnetic moments with 2% accuracy is obtained for m/7 = 1.818. For this 

value proton radius is more than experimental data Y(r^ = 0.870 ± 0.008 fm [9] and is 

equal to ^ {r\) = 1.291 ± 0.003 fm. For simultaneous description of experimantal data 
it is need to include negative contribution into the proton mean square radius, which is 
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equal to A(r|) = -0.910 fm^ for experimental data y (rl) = 0.870 ±0.008 fm. Assuming 
that this contribution is related to the valence quark radius, we have {r\^m^ = —1.875. 

Notice, that nonzero quark radius does not influence on the nucleon magnetic mo- 
ments, by using of which the ratio m/7 = 1.818 was obtained. If quark does not have size, 
then value, obtained in this model is corresponding to experimental data for m/7 = 0.4 
and m/7 = 0.9 (see figure a). For this m/7 it is possible to describe the nucleon magnetic 
moments but with less accuracy (m/7 — 0.4: deviations are 11% for proton and 20% 
for neutron), (m/7 = 0.9: deviations are 4.5% for proton and 11% for neutron); or in 
other case if we introduce large the quark anomalous magnetic moments, which should 
be different for u- and d-quarks [2]. So, if we assume that good description of the proton 
and neutron magnetic moments with the zero quark anomalous magnetic moments is not 
accidental, then model indicates that the valence quarks must have the negative radius. 

Conclusion. 

So, in this model the nucleon radius is defined by undimensional parameters m/7 and 
(rfg)m^. Absolute values of the quark characteristics (the quark mass and the quark 
radius) are not defined without specification of the "scale" -parameter, which defined the 
coupling energy. Oscillator wave function was used here for calculation, but the radius 
expression (16) with factor is general; it is possible that additional powers of m/7 

can change series (16) for other form of wave function and other definition of parame- 
ter 7. Dimensional factor 1/M^ is consequence of relativistic kinematic of the relative 
quark motion (see (8) and (9)). Scale invariant dependence of the nucleon radius and 
magnetic moments on m/7 is related to that the wave function and quasipotential are 
independent on bound system mass M in the impulse aproximation [1]. In the general 
case the relativistic expression for radius of the bound system contains infinite series of 
negative powers of m/7. In the oscillator model this gives large negative radius of the 
system with strong coupling of the spinor particles. Wigner rotation matrices gives small 
contributions into the proton radius in the region of parameters, which correspond to the 
the experimental data of the proton radius. In this model it is need to intoduce nega- 
tive quark radius for simultaneous description of the proton radius and nucleon magnetic 
moments. 
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